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Introduction . The earth's magnetic field is demonstrated 

i 

most simply by the forces exerted on a compass needle. The 
first man to appreciate the global character of the. geomag- 
netic field seems to have been William Gilbert. He deter- 
mined that the magnetic forces are guided from within tjie 
earth, and described his experiments in a treatise "De? TOagnete" 

published in 1600. • By observing the magnetic forces in the 

\ 

vicinity of a spherical lodestone and comparing them with what 
was known of the earth' s field, he concluded that the earth 
itself is a great magnet. A small bar magnet (magnetic dipole) 
of appropriate strength placed at the center of the earth 
would have the same external field as a uniformly magnetized 
sphere. Every schoolboy has mapped a dipole field by sprink- 
ling iron filings in the vicinity of a bar magnet. About 90 
per cent of the main geomagnetic field can be represented by 
a short earth-centered dipole inclined at an angle of about 
11.5° bo the axis of rotation. 

Although the compass was used in navigation during the 
fifteenth century, the earliest measurements of the deviation 
of the compass from true north were not made until the first 

half of the sixteenth century. The first world chart of the 

* ' 

variation of compass direction (declination) was constructed 
by Halley in 1701, with the knowledge that it would soon be 


obsolete 
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The earth's magnetic field is ever changing and its 
variations are recorded continuously at magnetic observ- 
atories over the globe. Nonmagnetic ships* aircraft* 
rockets and now earth satellites contribute their geomag- 
netic measurements to complete the description. Magnetic 
charts are currently revised every five or ten years tp ^ 
monitor the secular changes arising from fluid motions^* 
generating electric currents within the molten conducting 
core of the earth (see volume 5 for dynamo theory of the 
origin of earth's magnetic field). A World Magnetic Sur- 
vey has been undertaken during the years I964-I965 and will 
lead to improved world magnetic charts for the epoch 1965 * 0 . 

This project continues the highest traditions of interna- 
tional scientific cooperation established during the IGY. 
Measurements . In order to measure the vector geomagnetic 
field intensity it is necessary to determine a set of three 
independent quantities. Surface observatories usually meas- 
ure the horizontal intensity H, the vertical intensity Z and 
the magnetic declination D. With the introduction of extreme- 
ly accurate precession magnetometers to measure the magnitude 
of the total field B, the quantities determined are often B, 

D and the dip angle or inclination I.. It is sometimes con- 
venient. to use the set given by the north component X m H cos D, 

the east component Y m H sin D, and the vertical component 

.2 


H tan I* where B 


X 2 + Y 2 + Z 2 , 
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Figures 1 to 5 are reproductions of world magnetic 
charts of the surface field describing D, H, Z f B, and I 
respectively* The epoch and source are given for each iso- 
.magnetic chart. The isomagnetic lines of declination are 
called isogonic lines, and the zero isogonic line, along 
which the compass points to geographic north, is the agonic 
line. For the field intensity components the lines are 
called isodynaniic lines, and on the inclination chart the 
lines are isoclinic lines. The secular variation, which a- 
mounts to about 0.1 per cent of the total field per year, 
can be found on world isoporic maps prepared by the U.S. 

Navy Hydrographic Office in the 1700 series. 

Recent evaluations of the main geomagnetic field (Cain 
et al., 1965) have succeeded in reducing the errors in the 
field description to a fraction of 1% in regions of good 
data coverage. High speed computers are handling the data 
directly, thus eliminating the errors that may appear in the 
construction of intermediate charts. Surface anomalies on 
the order of a 100 km in extent do not appear on the charts, 
but at times may be a few per cent of the total field. 
Fortunately, the effects of these magnetic sources in the 
earth's crust disappear rapidly with altitude, and present 
no problem in describing the geomagnetic field in space. 
Pinole Field . The first approximation of the geomagnetic 
field given by a uniformly magnetized sphere is equivalent 
-to a centered dipole aligned parallel to the magnetic axis* 
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In Figure 6 let the magnetic axis be NS where N is the north 
geomagnetic pole and S the south geomagnetic pole, the radius 
of the sphere is a, the magnetic colatitude is ©, then at a 
distance r from the earth's center the magnetic potential 
V= -Mcosd/^2 where M is the magnetic moment. The intensity 


B of the field is given by B= -grad V 


The inward radial 


component of the field is *= 2Mcos© and the component 

perpendicular t<? r, in the direction of decreasing © is 


M 




£ . Therefore B = (H^+Z 2 ) 2 *= |*"l + 3 cos 2 0 **f 


1 

2 


r^0 

and since at the surface of the sphere at the equator M " a J Bo, 
this yields B B Bo (f-) 3 £1+3 cos 2 ©] ^ 

and tan I *■ % m 2 odt © . The dip in a dipole field is seen to 
n 

be the same at all points along the radius OP. The geomagnetic 
equator is the dip isoline 1=0 for the centered dipole field. 
The observed curve I=© for the actual main field is the mag- 
netic equator. The polar equation of a dipole line of force 
is r — c sin 0, where c is the distance at which the line 
crosses the equatorial plane. Clearly a field line that is 
at the surface at 0 B 30° crosses the equatorial plane at 4. 
earth radii from the center (see Figure 6). 

Spherical harmonic analysis . The first spherical harmonic 
analysis of the main field was made by Gauss in I838. Several - 


investigators using improved data and varying numbers of co- 
efficients have continued the work. The early studies are 
discussed, in Chapman and Bartels ( 1940 )• References to the 
recent work of Vestine, Finch and Leaton, and Cain may be 
found in Cain et al. (1965)* 
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With the assumptions of internal origin of the main 

magnetic field and negligible electric currents within the 

2 

magnetosphere, Laplace's equation V V=0 is satisfied and 
the magnetic field is derivable from a magnetic potential 
of the f.orm 

Co 

v - a 

n=l 

\ 

where 9 and are here the geographic colatitude and east 
longitude, Pjj(cosd) are the Schmidt normalized associated 
Legendre polynomials of degree n and order m and g™ n d 
h n are the Gauss coefficients determined from the surface 
magnetic data. The potential function is not a measured 
quantity, but the gradient of the potential function gives 
the magnetic field components that are observed at the 
earth's surface: 

X “ r 1 Z = 

rsin9 



t r*) n+1 R 

m^O V / n 


cos m 


X + h 


m 


n 


sxn 


3 p ” 


(cos9) 


In recent years, the most widely used coefficients have 
been those derived by Finch and Leaton, however the coeffi- 
cients of Cain and collaborators, determined by fitting the 
measured data rather than the data Tead from magnetic charts, 
are likely to be used more extensively in the future. 

The first approximation of the earth's magnetic field 
as a dipole can be represented mathematically by taking terms 

m 

in the potential function with ri=l (3 Gauss coefficients). 

This is the earth-centered dipole inclined at an angle of 11.5° 

/ 
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to the axis of rotation. The axis of the dipole intersect- 
ing the surface of the earth defines the geomagnetic poles 
and consequently the geomagnetic coordinate system. Figure 7 
indicates the geomagnetic coordinate system superimposed on a 
Mercator projection in geographic coordinates after Vestine. 

The geomagnetic poles are not the magnetic poles indicated 
on the surface magnetic field maps. The latter are the dip 
poles, which ai*e points on the surface of the earth where the 
magnetic field is vertical. The dip poles are not antipodal. 

A second approximation to the observed geomagnetic field 
is provided by taking terms through h=2 (8 Gauss coefficients). 

This is the eccentric dipole approximation and is provided by 
a suitable magnetic dipole parallel to the original central 
dipole, but displaced from the center of the earth about 
43(6 km towards a point at latitude 15.6° N, longitude 150.9° E 
for epoch 1955* The axis of the eccentric dipole intersects 
the eartli' s surface at points which represent another pair of 
magnetic poles. The approximate geographic location of the 

i 

various magnetic poles discussed above are shown in Table 1, 
with source and epoqh given . 

Table 1 

Geographic Latitude and Longitude of Various Magnetic Poles 

Observed Field 75.5°N, 259*5°E Epoch 1965 / 

(Dip Poles) 66 . 5°S, 139.9°^ " .Watford et al.(l965) ^ 

Centered dipole 78.5°N, 291. 0°E Epoch 1955 

(geomagnetic poles) 78. 5°S, 111.0 E Finch and Leaton(l957) 

eccentric dipole . 8l.0°N, 275.J°E ' Epoch 1955 

75*0 S, I20.4°fi Parkinson €T*Clea. y(19 58 ) 


More accurate approximations to the observed field can 
be made by taking additional terms in the potential expansion. 
The analysis by Cain et al. for the epoch 196 5 in Figures 4 
and 5 takes potential terms through n = 7(63 Gauss coefficients). 
Main field above the earth . At least 99 per cent of the main 
magnetic field on the surface of the earth originates from 
sources within the earth. When the magnetic field is known 
at all points on a closed surface which encompasses the sources 
of the field, the field is completely determined at all points 
outside the surface. The potential of the main field given 
above may be expressed in the form 

60 n+2 

V - Z V-r T " ' 

n-1 r 

where 

n m 

Tn = IET (g™ cos m?v + h m sin m^X ) P (cos©) 
m— 0 n n n 


m m 

and g and h„ determine the surface field. In order 
n n 

the field at height h above the earth, let a^ = a+hj 
be written in the convenient form 



to compute 
then V may 


r 

It is seen that the Gauss coefficients of degree n are dimin- 

n+2 

at height h. Local anomalies which 
correspond to harmonics of a considerably higher degree than 7 
are thus reduced effectively with elevation. The intensity of 


ished in the ratio' 


iof; 
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the dipole field and its components decreases outward as 
the cube of the distance from the center of the earth. The 
contributions to the geomagnetic field of harmonic terms of 
n“2, 3, 4, 5. ...decrease outward proportional to r” 4 , r 

— 6 _7 

r , r , respectively. Consequently, the greater 

the distance from the earth, the more like a dipole is the 
permanent field; provided that current sheets and ring 

currents are excluded from consideration. It may be expect- 

"> 

ed that the fit between computed and observed values of the 
permanent field will be better at altitude than at the sur- 
face, since the use of a limited number of harmonics cannot 
represent accurately the complex character of the anomalies 
in the surface field. 

Magnetosphere . For many years the region far above the 
earth’s atmosphere was thought to be a vacuum and conse- 
quently the earth's magnetic field was thought to extend 
outward in a completely predictable fashion. However, with 
the observations of comet tails by Biermann, theoretical 
studies by Chapman, and the magnetohydrodynamic theory of the 
solar plasma flux by Parker (see Solar Wind, volume l) , the 
permanent presence of streaming corpuscular radiation in 
space was established and verified by satellite observations. 
The conducting solar plasma confines the geomagnetic field 
within a cavity known as the magnetosphere, having a boundary 
called the magnetopause (see Magnetosphere, volume l) . Since 
the solar plasma is highly supersonic in the magnetohydro- 
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dynamic sense (velocity higher than the Alfven velocity in 
the medium) , a detached shock wave is produced in the medium 
ahead of the boundary. Viscous forces acting along the 
sides of the magnetopause are likely to be responsible for 
pulling out the magnetic field lines into>an extended geo- 
magnetic tail on the night side. (See Figure 3 of Geomag- 
netic Disturbances, volume 1 for view of magnetosphere) . 

Quiet day geomagnetic v ar iations . Examination of magneto- 
grams during geomagnetic ally quiet periods will readily re- 
veal characteristic magnetic variations that proceed accord- 
ing to local solar time and have amplitudes of approximately 
20 to SO gammas (gamma = 10 3 gauss) . This is the solar 
quiet day variation Sq. However, careful statistical analysis 
of the record is necessary to reveal another much weaker varia- 
tion (L) related to the lunar day and having an amplitude about 
one order of magnitude less than Sq . Analysis of both Sq and L 
reveal that about two-thirds of the magnetic variation is due 
to external sources, and one-third to sources within the earth. 
The internal contribution is likely produced by electric cur- 
rents induced in the earth by the varying external magnetic 
f ield . 

The quiet day solar variation of the horizontal intensity 
has an amplitude that is greater during the day than at night, 
and changes phase at a latitude somewhat below 4 O 0 . The am- 
plitude of Sq is greater in summer than in winter and is 
greater in years of sunspot maximum than in years of sunspot 


minimum. The lunar daily magnetic variation L also has a 
greater amplitude during daylight hours and in summer. The 
summer to winter ratio of L (2.6 to l) is much greater than 
for Sq . However, L increases only slightly with sunspot 
number, while Sq increases by about a factor of two from 
years of small to large sunspot number. 

The nature of the ionospheric current system that might 
produce the observed Sq variation or the observed L varia- 
tion can be computed independently of any theory indicating 
the precise mechanism by which the current system is pro- 
duced. Illustrations of such current systems may be seen in 
Chapman and Bartels (1940)* 

Dynamo Theory-Ionosphere . The most widely held theory to ex- 
plain the origin of the ionospheric current systems that 
would produce the Sq and L variations is the dynamo theory 
suggested in 1882 by Balfour Stewart. Tidal oscillations 
in the E region of the ionosphere move conducting air across 
the geomagnetic field inducing currents by dynamo action. 

The conducting air is the armature of the atmospheric dynamo 
and the geomagnetic field is the permanent field magnet in 
the analogy. The observed magnetic variations Sq and L are 
supposed to be caused by the magnetic fields of these iono- 
spheric currents (often called dynamo currents) . The analogy 
is far from exact, since in the presence of a magnetic field 
the conductivity of a partially ionized gas is anistropic. 

The mathematical development of the dynamo theory is given 


in Chapter 23 of Chapman and Bartels (1940) with recent de- 
velopments and references in a review by Fejer ( 1964 ). 

Two lines of investigation contribute supporting evi- 
dence for the dynamo theory. In one approach the wind vel- 
ocities and the conductivities are assumed and the resulting 
current systems with their magnetic variations are calculated 
In the second approach the observed magnetic variations are 
taken as a starting point and with conductivity assumptions; 
the electric current systems and then the wind velocities are 
calculated. The internal consistency of these calculations 
lend broad support to the dynamo theory. 

Even during magnetically quiet days the ever-present 
solar wind interacting with the magnetosphere will increase 
the magnetic field on the day side relative to the night side 
The rotation of the earth under this pattern will provide the 
observer on the ground with a quiet day magnetic variation 
related to currents flowing at the boundary of the magneto- 
sphere. Investigations along these lines are almost non- 
existent, and current explanations for the quiet day magnetic 
variation look toward the ionosphere and the dynamo theory. 
Some relationships of the geomagnetic field . l) In associ- 
ation with activity on the sun, enhanced solar* plasma inter- 
acting with the earth's magnetosphere causes geomagnetic 
disturbance and a host of related phenomena (see Geomagnetic 
Disturbances, volume l) . 2) The earth's magnetic field acts 
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on cosmic rays like a huge magnetic spectrometer. The Various 
particles are sorted with low energy particles reaching the 
polar regions* while only high energy particles are detected 
near the equator (see Cosmic Radiation, volume l). 3) Within 

the earth’s magnetic field, trapped electrons and protons 
spiral along the lines of force and drift around the earth to 
populate the radiation belts (see Van Allen . Radiation Belts, 
volume l) . ^ • 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


. World magnetic chart of the surface field: 
declination (D), epoch I960 

(Geomagnetism and Aeronomy, vol • 4 ( 1964 ), p.602) 


. World magnetic chart of the surface field: 
horizontal intensity (h), epoch I960 
(Geomagnetism and Aeronomy, vol. 4 ( 1964 ), p. 603 ) 

. World magnetic chart of the surface field: 
vertical intensity (Z), epoch i 960 
(Geomagnetism and Aeronomy, vol. 4 ( 1964 )# p. 604 ) 

. World magnetic chart of the surface field: 

total intensity (B), epoch 1965# after Cain et al. 
(JGR, vol. 70 (1965), p.3659) 

. * • 

. World magnetic chart of the surface field: 

inclination (X), epoch 1965 , . after Cain et al'. 

(JGR, vol. 70 (1965), p.366l) 

. Dipole Field 

* 

. Geomagnetic coordinates (curved lines) 

superimposed on a Mercator projector of the 
•world. 

(Handbook of Geophysics, Revised Edition I960, p.10-9) 
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data used were 250 determinations in Argentina during the period 1951-1954 [21, 22], 70 determina- 
tions in Colombia in 1958 and 60 determinations in Mexico [23]. - 

Asia. Recent observations of a high accuracy were used in Japan and the near-lying islands 
for 1951-1955 at more than 1000 stations [24, 23 J , in India for 1957-1958 and in the Philippines for 
1954-1955 [29] at 100 stations. We had no new data for the remaining territory of Asia. 

Africa. New determinations at 650 stations were used in North Africa [30, 31], in the Sahara 
[32-36] and in the Madagascar area [40] for 1947-1957. 

Australia. We used determinations made from 1927 through 1958 at 775 stations [4l]. Their 
distribution over the territory of Australia is uneven. Its eastern part is covered hotter. These 
data supplemented information on the magnetic field of Australia considerably. In the case of New 
Zealand we used 237 determinations for the period 1941-1951 [43]. 

Oceans. The magnetic determinations of the scientific expeditions aboard the schooner Zaryi 
were a considerable contribution to the study of the magnetic field of the oceans. These expeditions 
during the period from 1957 through 1961 made a survey in the Atlantic, Indian and Pacific Oceans, 
covering them with a thin but reliable network of magnetic determinations. The total length of the 
track followed by the Zarya was approximately 160, 000 nautical miles. The accuracy of determina- 
tions of declination w.as ~ 0 °. 5 ; accuracy of the horizontal and vortical components was ~ 0.001 oe. 

About 250 determinations of recent years [45, 46] were used for the islands situated in the 
Atlantic, Indian and Pacific Oceans. 

Arctio. In tho Arctic in recent years magnetic determinations have been made by expeditions 
on drifting stations and other high- latitude expeditiono. We used the results of these determinations 
for 1955-1960. (Earlier observations on drifting stations were used when compiling maps of the 
epoch 1955.) These determinations covered a considerable part of the north polar cone and were tbs 
basic data for map compilation. In addition, obsorvatlonal data for Spitsbergen for 1951-1953 and 
earlier were used [47 j . We had no other data auaea on foreiga observations and therefore the 
western sector of the Arotio was mapped solely on tho basis of old very modest and inexact data. 

$02 
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Fig. 2 


i Antarctica. This region, virtually unstudied in magnetic respects, was supplemented by 
fapietlc observations of very great importance for compilation of maps. These observations were 
•tdc by Soviet antarctic expeditions on sea and land; in the Indian Ocean determinations were made 
tirard the diesel- electric ship Lena in 1056. In the Pacific Ocean magnetic observations were made 
. &«rd tho diesel- electric ship Ob^ in 1957-1958 [48, 49], These observations were made along the 
1 ticka of the vessels to Antarctica. Their accuracy was low, but they were useful for such mag- 
1 Mtlcilly unstudied regions. 

I Soviet investigators made determinations along the coast of Antarctica and in the interior of 
I 4* continent. Although the number of determinations in Antarctica was not great each reconnais- 

i aace and each determination for that continent is of great value. The first traverse into the ant- 
irctic interior was made by P. K. Sen'ko in 1956 [50 ] . He followed a route about 100 km in length 
| fnm Mirny station to Pionerskaya station. In 1957 N.D. Medvedev made a traverse 1,400 Km in 
f injth from tho Pravda Coast to tho earth's South Geomagnetic Polo. Then, in 1958, determinations 
i *»ra made along the traverso Pionerskaya - Vostok I - Komsomolskaya - Sovetskaya - Pole of 
| WUtlvo Inacccssability [51} . In 1959-1960 determinations were made along tho traverse Vostok- 
? Klh Goographic Pole, and in 1960-1961 along the traverse Mirny - Komsomolskaya - Vostok and 
| atbo area of tho Shackloton and Wost Ico Shelvos. In addition, certain observational data of foreign 
f npedltlons wore used. Among thoso wore a survey of 1957 in the area of Prinoo Harald Land [52j, 
j •Uervatlonu of Amorican expeditions of 1951-1956 in tho area of the Palmer Peninsula and at 
* Lal« America and Little Rockford stations, French observations in Adelie Land [53] and at the 
f »hUonB Charcot and Dumont D'Urvtiie In 1951-1952 and Australian observations of 1953-1954 at 
k fcivion station and to the south of it [54] , 

[ All the data used'in compilation of the maps were reduced to the single epoch 1960. The 
| (vdoctlon was accomplished using world maps of isopora of D, H and Z for the period 1954-1959 
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[60 ], compiled by the Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation 
The cartographic base for compilation of maps for foreign areas was a grid at a scale of 
1:10,000,000. The isolines were drawn by interpolation between the value of the magnetic eiemt. 
at observation stations. 

For the area of the USSR we used the above-mentioned magnetic mapB, first subjected to k: 
simplification. o 

The principal interval between isolines was 0.005 oe for H, 0.01 oe for Z and 1 for D. Ai 
a result of the nonuniform magnetic study of the earth it is extremely difficult to make an evaluate 
of the world maps, expecially for areas which include the European part of the USSR, western ; 
Europe, the United States, Japan, India, Australia and New Zealand, it will be ~ 0.002 oe for H ■ 
and Z and 0^> for D. 1 

We have published maps obtained by replotting isolines from maps at a scale of 1:10,000,0# 
with the interval between isolines reduced to 0.01 oe for H and 0.002 oe for Z (Figs. 1-3). Tbo 
world magnetic maps for the epoch 1960 were prepared using more extensive and up-to-date dati 
than the maps for the epoch 1955 and therefore now are the most reliable magnetic maps of the wof 
On the oceans a great influence on the increase of quality of tho maps was the result of re c« 
observations on the schooner Zarya , whose contribution was great with respect to both quantity X 
quality. On ail the world maps compiled earlier the magnetic field for ocean areas was repress*' 
th»r mesh carh the hmhs nf -jhjfemtur.ca marie ca the r&vhi t/t fee 

ztw rresrr. The .x ~i Se.se fc the W-i 'A rv«e wap v/s* 


iia armirs ao Sa nsnlta and lowered aw accuracy of t*c 

nr tfa a™ C.W. m ** .* &« t»» ** *M* ft*** 
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OXSKN 

•,j wo soo only slight di(Toronc('s 
jjoch I>eing five yours later and 
el present charts use. f>3 instead 
efficients. The reasons for this 
'binge changes in the field during 
;jd are small and that the nmpli- 
ijtional l. r > coefficients are small 
fje of the lower terms. 

(smooth representation, it is of 
Ijle to illustrate the many small- 
fhat are found near the earth’s 
litours are thus niueh smoother 
^ standard world magnetic charts 
'illy constructed by other meth- 
rity of standard charts, however, 
■jmpt to show as much detail as 
with the result that there are 
S around anomalies in isolated 
d data coverage, but smooth 
ere. It is likely that such maps 
ire 3 give a more realistic csti- 
I averaged over the whole of the 

feature of the secular change 
j) in Tabic 4 is the fact that, the 
jot dominate the series as does 
iven though there are fewer 
terms determinable from the 
ing seeular change maps arc at 
jx as those for the field itself, 
ngo in two of the components, 
[n in Tigurc 4. The characteristic 
map is that the earth is almost 
: Ixitween positive and negative 
is pattern of positive Z is equally 
centers of negative change occur 
' and the positives are centered 
jitudes and over Asia. Since the 
hnces in the major secHdar change 
nailer than those for the cocffi- 
•in field (cf. Tabic 5), the details 
illustrated are equally less likely 

other models. To evaluate 
oh the presently derived set of 
res with other models, wc have 
, ations between the available 
liiputcd field components at the 
ation. Complete validation can 
fjincc it is not possible to make 
regions that arc void of data. 
Ijjta arc accrued it should always 






